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Abstract 

We present and study the results for the standard model process e + e~ — > vvbb at cm 
energies 150< yfs{GeV) < 240 and for Higgs boson masses 80GeU < m H < 120GeV, obtained 
from all tree-level diagrams and including the most important radiative corrections. The matrix 
elements have been calculated by the 'spinor bracket' method without neglecting masses, while 
the phase space integrals by an importance sampling Monte Carlo numerical integrator. The 
\/s dependence and the interference properties of the Higgs boson contribution and of various 
coherent background contributions to the total cross section are examined. The effect of the 
QED initial state radiative corrections is estimated. The important differential distributions 
for the Higgs boson and the background components are studied, providing information usefull 
for choosing cuts in Higgs searches. We also examine the effect of a minimal set of cuts and 
evaluate the importance of the WW fusion for detecting a higher mass Higgs boson at LEPII. 
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1 Introduction 



Increasingly precise LEP I data []]], combined with an evidence for the top quark of mass 
176 ± 15 GeV |§, have made feasible a very important —and impressively successful— test of 
the Standard Model (SM) to an accuracy level of one-loop quantum corrections 0,0. These 
data are also compatible with the Minimal Supersymmetric extension of the SM (MSSM |3j], 
which (at a price of introducting many more particles) provides perturbative stability against 
radiative corrections in a grand unification scheme. 

Yet, comparison to presently available data and requirements of theoretical consistency do 
not restrict sharply the mass of the SM Higgs boson, and even more so the parameters of the 
Higgs sector of the MSSM ||. On the experimental side, LEP I collaborations converge slowly 
towards a combined lower experimental bound m# > 65 GeV for the mass of the SM Higgs 
boson [|1[] : given this bound, the other results of LEP I are well-known to be rather insensitive 
on m/f. As for the theoretical constraints, the radiatively corrected SM Higgs potential V[<f)} 
develops an instability ||, and depending on the assumptions about the cure of this quantum 
catastrophe, some limits on m# are implied. E.g. assuming that the instability is avoided 
perturbatively by onset of new physics at A « 10 3 GeV (respectively A « 10 4 GeV) implies, 
for m t = 174 GeV and a s {m z ) = 0.118, a lower limit m H > 75 GeV (respectively m H > 100 
GeV) 0. A change Am t of the top quark mass induces a change Aran « 2Am t of such a 
limit 0. On the other hand, assuming the MSSM with m t ~ 170 GeV and supersymmetry 
breaking scale m s rs 10 3 GeV implies an upper limit < 125 GeV on the lightest mass Higgs 
boson, the limit changing with m t according to Aitlh ~ Am t M. The same upper bound (for 
m t ~ 170 GeV) was found from a low energy renormalization group analysis of the Higgs sector 
of the 'next-to-minimal' supersymmetric SM ||. Actually, because of the proximity of a fixed 
point in the low energy evolution of m t , one would expect that m# is well below this upper 
bound, expecting e.g. m# < 105 GeV for the above values of m t and m s [TU . 



Detection of at least the lightest Higgs boson (HB) and experimental study of its properties 
is presently among the highest priority goals in particle physics. Yet, LEP II will provide the 
only opportunity in the near future to search directly for a HB of mass m# > 65 GeV. It 
is thus very important to have complete and carefully analysed results for all e + e~-collision 
SM and MSSM processes with an 'uncoverable' HB signal at 170-210 GeV cm. energies. For 
the first time many of these processes will have 4 fermion final states, arising partly through 
2 intermediate heavy bosons, and we should learn to optimize the search for the HB in such 
processes and to push the explorable m H range as high as possible. A precise knowledge of the 
experimental background to these processes is, of course, also much needed. 

In this work we assume the validity of the SM at LEP II energies. Since the SM Higgs 
boson detectable at LEP II decays predominantly to bb, the most promising SM processes are 



e + e — > £ + £ bb, uubb, qqbb 
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where i = e or //, summation over the neutrino species is understood, and qq is a quark- 
antiquark pair generating jets just like the bb pair does. In the present study we concentrate 
on the complete tree level calculation, including approximatively the most important radiative 
corrections, of the process 

e + e~ — > uubb (1) 

and examine some problems of extracting the HB signal from it at LEP II despite a large 
background. We should mention that process (1) is of interest at LEP II for another reason: it 
will constitute an important background in the search for the process e + e~ — > xxbb, where x 
would be the lightest neutralino, stable when R-Parity is respected. Both z/'s and x's (if they 
exist) will be 'felt' only through energy-momentum conservation. 

An intermediate SM Higgs boson contributes to process (1) through the so-called Hig- 
gsstrahlung and WW fusion diagrams of Fig. 1. This HB signal is embedded in a coherent 
electroweak background arising from the 20 diagrams fo Fig. 2. These can be classified as 

• the e-exchange diagrams of Fig. 2a, where ZZ or Z7 coexist and cascade on the s-channel 

• the WW multiperipheral and fusion diagrams of Fig. 2b 

• the W- and e-exchange diagrams of Fig. 2c, with a radiation of 7 or Z 

• the s-channel e + e~ annihilation diagrams of Fig. 2d, with Z radiation from a final fermion. 

For reasons of brevity, the background processes of Fig. 2a, b, c, d will be called respectively 
ZZ/^j cascade, WW fusion, W exchange and Z radiation processes. We work in the t' Hooft- 
Feynman gauge, but have ignored diagrams with a coupling of the Higgs boson or of Goldstone 
bosons to e ± , u e , V e leptons, since they contribute insignificantly. 

Experimentally, one must extract the HB signal form all final states consisting of a pair 
of imperfectly 6-tagged jets + missing products. The important components of the incoherent 
background, denoted Bq CD , Bww, Bw, Bee and B^ for brevity, correspond to the following 
processes: 

• Bqcd'- s-channel (virtual 7/Z mediated) production of bb or cc jets (there is the possibility 
of misidentification of c jets in 6-tagging), accompanied by soft undetected hadrons or by 
beam radiated photons escaping detection. 

• B ww : WW production, with one W decaying to cs jets and the other decaying undetected 
through a chain W — > vr — > uV+ soft charged particles. 

• Bw' Wev events, with e escaping undetected and W decaying into cs jets. 

• Bee, t = e or ji: £ + £~bb, £ + £~cc final states, with the leptons escaping in the forward 
and/or backward direction. 
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It is evident that very efficient criteria of event acceptance must be invoked in order to 
isolate (if possible) the Higgs component at a detectable level. 

In order to put the present study in a proper context, we recall briefly previous work 
on e + e~ — > vvbb. Starting from , there have been several ambitious simulations of this 
process by experimental groups, aiming at a realistic estimate of the SM Higgs boson signal 
to background ratio at LEP II (e.g. Hl2| , fl3| ) . Yet, they ignored the WW fusion diagram 
assuming that it was important only at higher energies. But as shown by our 

(a) it increases 



of Fig. 1, 

results, WW fusion is already important at LEPII energies for two reasons: 
the signal by about 10% around its maximum and (b) it increases by large factors the weak 
signal around the kinematical threshold for HZ production. They also ignored parts of the 
electroweak background and evaluated the rest approximately, invoking only 2^2 body 
processes followed by appropriate decays. Theorists, on the other hand, were contended to 
study e + e~ — > vvH without quantitative studies of the background. They usually included the 
WW fusion contribution to this 2 — > 3 body process |T4] , [T5| , [|TBJ , but stressed its importance 
only for ^/i > 500 GeV [17|,[j8]- It was only recently |]19|,[|2(| when our present work was well 
under way, that the process e + e~ — > uubb was calculated including all signal and background 
contributions, with emphasis on results for -y/s above 200 GeV. 



Despite some overlap between |L9]],|2(J and the present study, we decided to complete and 



publish our results for several reasons. We shall not dwell on a trivial one, the need to check in- 
dependently calculations involving 22 diagrams and multidimensional phase space integrations, 
since we find good agreement when calculating same things with same parameters. More impor- 
tantly, in view of the necessity of event selection criteria for improving the signal/background 
ratio, we put emphasis on studying the detailed characteristics of the signal and the coherent 
background, including some differential cross sections not considered in [O,p0f. This informa- 



tion can be beneficially combined with experience about Bq CD , B ww , B w and B u to choose a 
set of experimental cuts aiming to extract a measurable HB signal. We also introduce several 
quantitatively important technical improvements over previous work (approximative inclusion 
of the QED initial state radiative (ISR) corrections and of important QCD radiative correc- 
tions to the adopted HB width and the Higgs— bb coupling, approximative employment of an 
Improved Born Approximation), and examine the effect of a set of minimal cuts. These results 
allow estimates of the SM Higgs boson detectability at LEP II more realistic than previously. 
As a by-product, we reevaluate the importance of the W — W fusion for detecting a higher mass 
HB at LEP II, already noticed in 
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f20|j . We should mention once more that, for reasons of 
experimental accessibility, all cross section presented below are summed over neutrino species, 
(in |2{| cross sections with electronic neutrinos are emphasized). 
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2 II. Method of calculation and computation 



The process e + e — > vvbb interrelates six external spin ~ fermions. For fixed momenta of these 
fermions, the initial-spin averaged and final-spin summed matrix element squared is given by 



F 



(2) 



where A runs over a basis of spin states for external fermions and F runs over all Feynman 
diagrams contributing to the process. Being interested in the Higgs signal, we must avoid 
approximating the b quark as massless: so we decided not to use this approximation for the 
electrons either. We then choose a convenient spinor basis in accordance with pll,[p2| and 



calculate analytic expressions for all iM x : the remaining operations in (2) are performed, for 
each point of phase space, by a computer program. (The advantages of this procedure in 
comparison to the traditional 'trace technique' are discussed in [^TJ,|22[). Before describing 
the structural form of iM[, we would like to briefly motivate (in a novel way) and present the 
spinor basis used, since it deserves —due to its usefulness— to be more widely understood and 
used. Its basic advantage is to extend smoothly to massive fermions the natural spin basis for 
massless fermions. 

A well-known covariant way of characterizing spin states is by invoking the translation 
—invariant Pauli-Lubanski axial 4- vector operator p3[. 



W fi P IM = 0, [W^,P u }=0 (3) 
[W», W v ] = ie^ xp P x W p , [W 2 , W v \ = (4) 

where are the spin-parts of Lorentz generators, W 2 = W^W^, and we use the conventions 
£ oi23 _ j anc j _ _ ^ _ ^ _ _ ^ gi ves j- p ^ w hen acting on e T 1 Pfl x M , i.e. on 

particle (antiparticle) solutions of 4-momentum P^. For P M timelike (massive particle), is 

spacelike and reduces in the particle rest-frame to the spin 3-vector E 

W = (0, ±mS) for P M = (m, 0) 

while for P^ lightlike (massless particle) W p is lightlike 'parallel' to P M , W p = ±/iP M , where the 
helicity h = P ■ W/(P ) 2 is a Lorentz-invariant characteristic of a massless particle. According 
to (3), (4) we can choose simultaneous eigenstates of P^, W 2 and of the projection of along 
some (eventually P^-dependent) 4-vector a spin basis is characterized by the choice of S p . 
This choice is best understood if we complement P^ 1 by three independent 4-vectors so as to 
form a complete basis of Lorentz 4-vectors convenient for expressing W* 1 . 

A P^ either lightlike or timelike can be complemented by two unit spacelike 4-vectors e^, 
such that they form an orthogonal sub-basis (a ■ b and a 2 mean a M 6 M and a p a p respectively) 

P • e ( i) = P • e (2 ) = e ( i) • e (2) = 0, e 2 {1) = e\ 2) = -1. (5) 
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The choice of the fourth basic 4-vector useful for both P 2 = and P 2 > is motivated as 
follows. 

Since for massless particles the helicity h (identical for such particles to the chirality) is 
a Poincare invariant property conserved by vector and axial-vector interactions, we choose 
eigenstates of h as a spinor basis for massless fermions. Thus we choose = P M for such 
particles. Complementing now P M , e^, efL by a lightlike 4-vector k^ such that 

k ■ e ( i) = k ■ e (2) = k 2 = 0, P-k^0 (6) 

we can write 

W ■ k 

W = ±hP^ = -—-rP* 1 (7) 
a form convenient for generalization to massive particles. 

Turning now to massive particles, P M becomes timelike and W 1 spacelike (W 2 = — m 2 E 2 ). 
Expressing such a in our Lorentz basis we find (using W ■ P = and (5), (6)) 



W ■ k 



rn 
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Thus choosing for massive particles 

TP 2 

S» = p» - (9) 

we obtain a spin-characterization which connects smoothly the massive to the massless case: 
as p 2 = m 2 _> o, W -> ±/iP", S» -> P^ and (8) -> (7). This ^ satisfies 

,S 2 = -m 2 , S-P = 0. (10) 

If /c M transforms to £^ when Lorentz-boosting to the rest frame of the particle, S 11 transforms 
simultaneously to (0, — ml), where I is the unit 3- vector along I, while transforms to 
(0,±mS). Hence 

P-k m 2 K J 

i.e. our spinor basis for a massive particle corresponds to eigenstates of the projection of the 
spin £ along the unit vector I in the rest frame of the particle. 

The final result of this procedure is to use, for both massless and massive fermions, as spinor 
W ■ k 

basis the eigenstates of — , where k^ is some fixed lightlike 4-vector such that P ■ k ^ 0, and 

P • k 

P M , are the 4-momentum and the Pauli-Lubanski 4-vector of each particle. It is useful to 
notice that for spin — fermions 

- \e^ xp Y, vX = i l5 E"> (12) 



(with the convention 75 = Z7°7 1 7 2 7 3 ), so that 

w - k = - m = l5 m -p-q (13) 

using YT V = |[7 M , 7 V ] and the notation = 0^7^ for any 4- vector a^. In this case our spin- 
characterizing operator can be written 

where + (— ) corresponds to positive (negative) frequency solutions, i.e. to particle (antiparticle) 
wavefunctions u(p) (v(p)) respectively. 

W ■ k 

Denoting the eigenvalues of — — — by A/2 (A = ±1), let U\(p, m), V\(p, m) be the correspond- 

ST ' rh 

ing fermion and antifermion eigenstates, solutions of the Dirac equation with 4-momentum p M . 
It is easy to see that, for both massive (m > 0) and massless (m = 0) particles, all the above 
requirements and the normalization-fixing relations 

U\(p, m)u\(p, m) — + m, v\(p, m)v\(p, m) = p 1 — m (15) 

A A 

are satisfied by the choice [^TJ 

u\(p, m) = -(j> + m)u-.\(h), v\(p, m) = U-\(p, -m) (16) 

T) 

where 77 = y/2p ■ k and 

$u x (k) =0, u x {k)V x {k) = \(l + \ 1 ^ (17) 

i.e. v~(k), v+(k) are the helicity eigenstates of a massless spin - fermion having the reference 

4-momentum k^. We can take e.g. v+{k) = fa-(k) where e M is any 4- vector satisfying e 2 = — 1, 
e ■ k = 0. One should notice that (17) corresponds to a normalization 

V x {k)-fu x {k) = 2P (18) 

which, in collaboration with the I/77 factor in (16), leads to 

u\(p,m)u\(p,m) = 2m, v x (p,m)v x (p,m) = —2m. (19) 

According to our spin specification 

u x (p,m)u x (p,m) = -[1 + A7 5 ^( 3 )](^ + m) 

= \ 1 + A ^(l"^) V + (20) 



where is the dimensionless unit spacelike 4- vector S^/m. The corresponding expression for 
v\(p,m)v\(p,m) is obtained by changing A — > —A, m — > — m in the r.h.s. of (20). Formulae 
(15), (16), (19) and (20) have a smooth limit for m — > 0, and give then the correct relations for 
massless fermions. 

Any tree Feynman diagram is characterized by a number of uninterrupted fermion lines, 
each streaming a fermion nb. flow which enters and leaves the diagram. Along such a line there 
is a contraction of Dirac indices of a sequence of spinors, interaction vertices and numerators 
of fermion propagators. Vector or scalar bosom propagators provide contractions of Lorentz 
indices of vertices belonging to different fermion lines. The essence of the 'spinor bracket' or 
'pseudohelicity method' of [ 21 1 , [22] is to express each such fermion line sequence in terms of 
few basic blocks. To exhibit the structure of each iM[ in terms of the building blocks of the 
present approach, substitute in the corresponding expression all numerators j> ± m of fermion 
propagators according to (15) and write all coupling constants as ClPl + CrPr in terms of the 

chiral projectors Pl = -(1 — 75), Pr= ~(1 + 75). Then each iMf becomes a sum of products, 
these last having as factors two types of complex scalar functions: the spinor brackets 

Wai(pi, rni)(C L P L + C R P R )ux 2 (p 2 ,m 2 ) 

(whose form motivates the choice of building units, see below), and the so called Z-functions 

Z(l,2,3,4) = [ux 1 (pi,m 1 )r i ux 2 (p 2 ,m 2 )][ux 3 (p3,m 3 )Tl l ux 4 (p4,rn 4 )] 

= r(c L p L + c r p r ), r; = ^(C L p L + c' R p R ) 

Decomposing next U\(p J ±m) into chiral spinors a\(p), v\{k) 

I 771 

u x {p, ±m) = -j)v_ x {k) ± — v-x{k) = a x {p) ± iw~\{k) 

7] 7] 

and using the so-called Chisholm identities for such spinors 

[aA(<?i)7 M "A'(<?2)]7M = 2$\\>{a\(q2)ax(qi) + a-. x {qi)®-\{<l2)) 
[ax(qh^\'(k)]'j tI = 2Sxx'Wx(k)a x (q) - a-\(q)V-x(k)] 
P\{k)l^ax'{q)}^ = 25 X \'[ax{q)Vx{k) - is_x(k)a_ x (q)} 



one can express straightforwardly, and once for ever ||21||,||22||, also the Z-functions in terms of 
the following basic units (here we slightly refine previous work) 

Yx ,R (pi,P2) =ux{pi)Pl,rUx{P2) 

Sx' R (Pl,P2) =Ux{pi)P L ,RU-x{P2) =U\{Pl)PL,Ra-\(j)2) 
TTl I 

\l = — , ce\(p)i/-\(k) = Vx{k)a-x{p) = J2p ■ k = 7] 

7] v 

The main properties and the explicit expressions (after adopting f+(k) = fcv-(k)) of the Y- 
and S'-functions are 

^a L (Pi>P2) = Y R x ( Pl ,p 2 ) = Y L x (p 2 , Pl ) = Yx R (p 2 ,Pi) 
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y L ipiiP2) = y R (pi,P2) = 

S+(pi,P2) = ^[(Pi ■ k )(P2 ■ e) - (pa • fc)(pi ■e)+ie(k,e,p 1 ,p 2 )] 

s R ( Pl , P 2) = -s^( Pl ,p 2 y, s*( Pl ,p 2 ) = s L ( Pl , P2 ) = 

where e(k,e,pi,p 2 ) = e ilv \ P k il e v p\p p 2 . l £0123 = — 1- We shall denote the nonzero S'-functions 
simply S±(j>i,p 2 ), neglecting the upper L,R index. For completeness of presentation, the 
expressions for Z-functions are tabulated below. For and e M we used the convenient choice 
k^ = (1, 1, 0, 0), e M = (0, 0, 1, 0) : k^ should be non-orthogonal to all external p^. 



Z(+, +, +, +) = 2[S + {p 1 ,p 3 )S^{p A ,p 2 )C R C' R + ^ 2 t]zT]aC l C' r + r]iV2^3^C R C' L ] 
Z(+, +, -, -) = 2[S + (p 1 ,p i )S-(p 3 ,p 2 )C R C' L + ihpl^tIaClC'l + r]iV2^3^C R C R ] 
Z{+, -, -, +) = 2[//i772A«3?74ClC^ + t]i^2V3^CrC'l ~ ^iV2V3^C L C' L - r]ifj, 2 fj l3 r]^C R C' R ) 

Z(+, -,+,-) =0 
Z{+, +, +, -) = 2r) 2 C R [n 3 S + {pi,p±)C' L + mS+fahPi)^ 
Z(+, +, -, +) = 27 ]l C R [fi 4 S4p 3 ,p 2 )C' L + li 3 S^(p 2 ,p 4 )C' R ] 
Z(+, -, +, +) = 2r] 4 C' R [fi 2 S + (p 1 ,p 3 )C R + fiiS + (p 3 ,p 2 )C L ] 
Z{+, -, -) = 2r] 3 C' L [fx 1 S + (p i ,p 2 )C L + fx 2 S + (p 1 ,p i )C R ] 

The remaining Z are obtained by exchanging simultaneously + <-> — and L «-> i? in the 
above expressions. 

To recapitulate, using the described spinor basis and following the outlined calculational 
procedure, each Feynman amplitude zMjf was expressed as a sum of products of the above 
defined Y-, S- and Z-functions. The Z-functions are expressed in turn in terms of the truly 
basic building units, the S— and F-functions, having simple properties and explicit expressions. 
The whole scheme is very convenient and allows compact programming. 

We use the t'Hooft-Feynman gauge and Breit-Wigner propagators for the Z,W and H bosons. 
The mass of W (m w = 80.23 GeV) and the mass (m z = (91.1888 GeV) and width (T z = 2.4974 
GeV) of Z have been taken as inputs; the width of W being calculated from the SM tree-level 
formula. 

The HB mass was varied in the range 80-120 GeV. For each m# the total is estimated 



using the program kindly provided by the authors of [p4|, which includes the significant QCD 
and electroweak radiative corrections to the tree-level width formula. The widths Th and F b i 
found for various m# are as follows: 

In order to include approximatively the corresponding QCD correction to the Hbb vertex 
in the Feynman amplitudes, we employ in the Hbb coupling a running mass mt)(m 2 H ) which, 
when used in the tree-level formula for T b i, reproduces the above stated widths. In the range 
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m H (GeV) 


80 


90 


100 


110 


T H (MeV) 


1.734 


1.911 


2.106 


2.393 


T bh (MeV) 


1.513 


1.661 


1.806 


1.948 



of masses considered here, rrib(m 2 H ) differs from the QCD running mass at Q 2 = m 2 H [JJ5] by 
a factor 1.12 ±0.01. The value of as{m 2 z ) was taken 0.123. The kinematical quark masses 
needed in various parts of our calculations are taken to be m;, = 4.7, m c = 1.45 and m t = 175 
GeV, while the KM matrix element Utb = 1- The gauge couplings were fixed by adopting 
a eff = 1/128.07, (1/137.0359895) at Z(j) vertices and using the effective g 2 = A^/2G F m 2 w 
{Gp—1. 16639565 10~ 5 ). These effective parameters partly include EW radiative corrections in 
our tree-level amplitudes, in the spirit of the 'Improved Born Approximation' [2C]. 

The QED ISR corrections were calculated approximatively using the universal radiation 
factor, incorporating virtual and soft contributions up to second order in a, with infrared 



photon exponentiation, and a hard bremsstrahlung contribution [27]. Our cross sections should 



not be sensitive to final state radiation, because vv do not radiate and bb are massive enough. 

Denoting final u, V, b,b as 3, 4, 5, 6 respectively, their Lorentz-invariant phase space volume 
element ciLips (s; p^PiP^Po) was decomposed into a product of effective 2-body elements 

dLips (s;p 34 p 56 )dLips (s 34 ; p 3 p 4 )dLips (s 56 ;p 5 p 6 )^- 

where 

fi jit , /U rp2 -Q. 

Pij — Pi < Pj > s ij — ^ij ~ Pij 

and each 2-body d Lips was evaluated in its cm. frame. Smoothing the integrand by well-known 
changes of variables |28|, the phase space integration was performed by the event generation 



and Monte Carlo integration program VEGAS [[^J, with estimated accuracy better than 1%. 

The fortran program exists also in a MC generator form, where it is interfaced to JETSET 
|34| and so after hadronization produces final state particles. [] The MC generator is controlled 



by a set of cards where one can determine the HB mass, the cm energy, the inclusion or not of 
radiative corrections, possible cuts on m v9 and m b i, the number of requested events, and whether 
one wishes to write the events in LUND format on some separate file. One can restrict, through 
switches in the same cards, the calculation to a subclass of diagrams and to the production of 
v e or only. Finally, all input constants used in the calculations: mz, raw, ^z, Gp, a (at 
and Mz), as and the quark kinematical masses, are given from the outside on the same set of 
cards, providing to the user complete control of the input variables for comparison purposes. 



The code is available from the authors upon request 
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3 Results and comments 



Considering SM processes e + e~ — > bb + ff at LEPII, it is well known that about 20% of 
the Higgs signal belongs to the final state vvbb, imbedded in a large background. Using the 
calculational procedure and the parameters presented in the previous section, we now proceed 
to investigate the characteristic properties of the HB-mediated (Fig. 1) and the coherent 
background (Fig. 2) contributions to e + e~ — > uVbb (sum over neutrino species), in order to 
provide usefull information for Higgs searches at LEPII. 

To gain an orientation we start with Fig. |||, which shows, for m# = 90 GeV and over 
cm. energies 150 < (GeV) < 240, the total expected number of e + e~ — > vvbb events 
for integrated luminosity 1000 pb -1 , together with the expected signal-only events from only 
the HB-mediated diagrams of Fig. 1 and the expected background events from the diagrams 
of Fig. 2. The expected signal events when one includes also the interference term between 
diagrams 1 and 2 are also plotted, but are indistinguishable (differences below 0.2 %)[] from 
the signal-only curve. In the same figure the two lower curves provide the signal-to-background 
ratio and the signal-to-Higgsstrahlung ratio, the second indicating the relative importance of 
the WW fusion contribution to the signal. It is interesting to note that the last contribution 
provides an important extension of the HB detectability limit at LEPII. We see from the figure 
that for ^fs clearly above the HZ threshold the WW fusion contributes about 10% of the signal, 
but for a/s around and below m# + mz its relative importance grows rapidly with decreasing 
cm energy. We will examine further down a realistic numerical example. 

Fig. [| shows the Higgs signal for m H = 80, 90, 100 and 110 GeV and the background 
cross sections at LEPII energies, before and after taking into account ISR corrections. These 
corrections depend strongly on y/s, ranging from 5% to 25% and being largest (as expected) at 
those energies where the cross sections increases most rapidly with ^Js. 

We next examine the relative importance of various background contributions and their 
energy dependence, paying particular attention to Eq M = 175, 192 and 205 GeV, which are 
the nominal energies studied in [Q. Fig. | shows the separate contributions due to the 



background processes of Figures 2a, 2b, 2c and 2d, without and with the ISR corrections 
respectively. The ZZ/ r y cascade contribution (Fig. 2a) dominates. It grows rapidly at the 
ZZ production threshold, while the other contributions increase more gradually with ^fs. At 
\fs = 175, 192 and 205 GeV the total coherent background cross section with ISR corrections is 
about 0.006 pb, 0.045 pb and 0.065 pb respectively. It is clear that the diagrams of Fig. 2d can 
be practically ignored at LEP II energies, and those of Fig. 2b contribute too weakly to provide 
some sensitivity on the exact value of m t . An interesting effect is a very strong destructive 
intereference between the WW fusion and W exchange diagrams of the background, as clearly 
demonstrated on the same figure. This is a reflection of the unitarity (through gauge symmetry) 
of the SM, which prevents these background components from violating asymptotically the 
log(s) bound. 

2 The smallness of this interference term was checked to hold for all Higgs masses including mn = mz 
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To decide about acceptance cuts, it is important to know the more accessible differential 
cross sections of the background and the signal contributions. As such we studied the differential 
cross sections versus the invisible invariant mass m^, the visible invariant mass m bb , the cos8 bb 
of the total visible 3-momentum p bb and the cos 9 b of the 6-quark 3-momentum in the collision 
cm. system. All differential cross sections have been calculated with ISR corrections. The 
above differential cross sections for the signal (m#=90 GeV) and the total background at 
x/s = 192 GeV are presented on figure || The next figure |7] shows how the background 
changes when decreasing y/s from 192 to 175 GeV. And Fig. |8| shows how the background 
differential cross sections split into approximatively gauge independent components ZZ/7, WW 
fusion +W exchange and Z radiation. 

The gross features of the background can be understood from the structure of the ZZ/j 
cascade diagrams, which dominates it. The sharp peak around m vv « mz is evidently due to all 
vV pairs being produced through an intermediate Z . The two peaks at m b i > and m bb ~ mz, 
are due to the 7 — > bb and Z — > bb cascades in the diagrams of Fig. 2a. The peak at m bb rs mz 
dominates as the energy increases and we pass the ZZ production threshold. One should also 
notice that the ZZj^i cascade and the W exchange contributions peak along the beam axis. 
This is due to the t- and w-channel exchange of a light fermion in diagrams of Fig. 2a and 2c, 
leading in the differential cross section to factors of the form (sin 2 6 bb + (mass) 2 /,?)" 1 , where 
(mass) 2 /s is very small for Z7 cascade and small for ZZ cascade (except for \fs just above the 
ZZ production threshold): for a discussion see e.g. p0 |. 



The peak of the signal around m vV « mz is due to the dominance of the Higgsstrahlung 
diagrams above threshold, while the WW fusion produces a broad m V y spectrum. The narrow 
width of the HB is responsible for the spectacular peak at m bb . One should also notice that, in 
contrast to the background, the signal angular distributions are remarkably isotropic. For the 
Higgsstahlung contribution this is in agreement with the analytic result |3l[ that the differential 



cross section of e + e — > HZ has angular distribution (in the cm. system) 

Q 2 sin 2 9 



2 



2m^ 

where Q = A 1 / 2 (s, m|, m 2 H )/2\fs is the 3-momentum of Z or H: for yfs just above m# + mz, 
Q 2 /2m 2 z is very small and the distribution is flat, but this will change sufficiently above the 
HZ production threshold. The angular flatness of the WW fusion signal is in accordance with 
the expectation of an isotropic distribution of the 3-momenta of the Higgs bosons produced by 
WW fusion at energies not far beyond the WW threshold. It is also clear that there should be 
considerable constructive interference between the two signal components for y/s around and 
below the HZ threshold. Indeed an examination shows that at these cm energies about half 
the increase of the signal cross section associated with the presence of WW fusion comes from 
the interference with the Higgsstrahlung, while at higher \fs this interference dies out. Similar 



conclusions have been reached in the very recent study of |32fl 



Knowing the main properties of the signal and the coherent background, one should next 
proceed to cope with the crucial problem of the very large incoherent background, whose main 
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components, as discussed in the introduction, are Bq CD , Bww, By/ and Bu (I = e and /i). This 
problem is clearly beyond the scope of this work: it has been studied by the LEP experimental 
groups as reported in We shall simply adopt the following event selection criteria, which 
were prescribed by the LEP200 event generator group, and study how they affect the signal 
and the coherent background: 

(i) 50 GeV< m bi (ii) (m z -25 GeV) < m uP < (m z +25 GeV) 

Experimental suppression of all background to the HB signal in e + e~ — > bb + missing 
products requires of course several further event acceptance criteria, such as allowing events 
with only two jets, both 6-tagged and of constrained morphology (acoplanarity and acollinearity 
cuts), forbidding evens with too small number of charged tracks or with large missing energy 



in the beam pipe, etc 33 



The total signal and background cross sections after cuts (i) and (ii), for various HB masses 
and without or with ISR corrections, are shown on Figs. [10|. Fig |] shows how these cuts 
affect the background components. We see that the coherent background is reduced greatly 
(by about 65%) at y/s = 175GeV, but the reduction is much smaller (10-15%) at y/s = 192 
and 205 GeV. This can be attributed mostly to cut (i) which supressed the photon mediated 
bb production, whose relative importance in the total background cross-section decreases with 
increasing y/s. As regards the signal, it is mostly influenced by cut (ii), which reduces the 
contribution induced by the WW fusion and characterised by a broad m vD spectrum. So cuts 
(i)+(ii) lower the signal by only 5-10% for y/s > rrin + Tnz, but the reduction grows rapidly as 
y/s decreases below the HZ threshold. 

It is interesting to estimate carefully, including the effect of our cuts, how the HB detectabil- 



ity limit at LEPII is extended by the WW fusion process. To this end, we show on fig 11] the 
y/s dependence of the (ISR corrected) signal cross section before and after cuts, for m H = 100 
GeV, and the corresponding total signal/Higgsstrahlung ratio. We see that for y/s clearly above 
rriH + rnz, the WW fusion contribution (about 10% of the signal) is reduced by about 30-40% 
by our cuts. In the region 180-200 GeV, where the WW fusion relative importance increases 
rapidly with decreasing y/s, our cuts reduce by only about 20% the signal-to-Higgsstrahlung 
ratio. E.g. at -^5=192 GeV about 50% of the events are due to WW fusion, this reducing to 
about 40% after the cuts. The final number of expected events for an integrated luminosity 
(e.g. by summing 4 experiments over 2 years) of 1000 pb" 1 is 8. 

Finally Fig. [L^ presents information which allows an immediate estimate of the reduction 
of the number of background and signal events due to the limited acceptance of micro-vertex 
detectors used to tag the b jets. It is clear that signal and coherent background are about 
equally lowered by this effect. Micro-vertex detectors extending to cos6 l =0.9 permit the loss 
of only about 10-15% of the events when b-tagging of both jets is required, while one has a 
negligible loss (around 2%) if contented with only one b-tagged jet. 

In order to allow comparisons with other existing Higgs codes, we followed the LEP200 
Higgs event generator prescriptions concerning the input values and obtained the cross-sections 
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requested for comparison purposes. We present them in the tables below. As can be seen in 
|33| the agreement is satisfactory, differences are of the order of 1%. 



m H (GeV) 


65 


90 


115 


oo 


175 GeV (i/^) 


63.447(0.289) 


2.316(0.011) 


1.257(0.006) 


1.257(0.006) 


192 GeV (v^) 


72.01(0.317) 


46.840(0.093) 


19.575(0.087) 


19.453(0.087) 


205 GeV 


67.704(0.301) 


53.569(0.240) 


25.531(0.114) 


23.646(0.105) 


175 GeV (u e u e ) 


70.103(0.302) 


4.996(0.023) 


1.051(0.006) 


1.051(0.006) 


192 GeV (u e u e ) 


79.158(0.348) 


52.744(0.140) 


20.605(0.102) 


19.678(0.098) 


205 GeV (v e v e ) 


76.165(0.359) 


61.429(0.280) 


31.048(0.146) 


25.737(0.122) 



Table 1: The process e + e~ — > vvbb at Born level, following the LEP200 generator Higgs 
group input values and prescriptions. Cross sections in fb.The numerical integration errors are 
included in parentheses. 



m H (GeV) 


65 


90 


115 


oo 


175 GeV (i/ M P„) 


54.186(0.283) 


1.660(0.008) 


0.917(0.006) 


0.917(0.006) 


192 GeV (^z/ M ) 


64.405(0.246) 


35.367(0.109) 


14.68(0.087) 


14.60(0.087) 


205 GeV (u^) 


63.826(0.222) 


45.736(0.171) 


21.275(0.105) 


20.058(0.105) 


175 GeV (u e D e ) 


60.806(0.294) 


3.663(0.018) 


0.775(0.006) 


0.775(0.006) 


192 GeV (u e u e ) 


71.151(0.268) 


40.212(0.164) 


15.222(0.098) 


14.61(0.098) 


205 GeV (v e v e ) 


71.794(0.266) 


52.543(0.199) 


25.428(0.124) 


21.76(0.122) 



Table 2: The process e + e~ — > vvbb including ISR corrections, following the LEP200 generator 
Higgs group input values and prescriptions. Cross sections in fb. The numerical integration 
errors are included in parentheses. 
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4 Summary and conclusions 



In this work we studied the SM process e + e~ — > uubb (sum over neutrino species) at cm 
energies y/s— 150-240 GeV and for Higgs boson masses m#=80-120 GeV. The results were 
obtained from all tree level diagrams and included approximatively the most important radiative 
corrections, i.e the QED ISR corrections, the QCD and EW corrections to the Higgs width and 
a QCD correction to the H66 coupling. The matrix elements were calculated by the 'spinor 
bracket' method without neglecting masses, and we tried to make the method more accessible 
by describing it in detail and motivating its spinor basis in a novel way. The phase space 
integrals were calculated by an importance sampling Monte Carlo numerical integrator. 

The signal is composed of the Higgsstrahlung and the WW fusion amplitudes. The first 
dominates for y/s above the HZ threshold, but the second becomes increasingly important as 
y/s decreases around and below that threshold, extending probably (depending on luminosity) 
to van > \/s — mz the Higgs detectability limit at LEPII. The two amplitudes behave similarly 
over phase space, providing in particular isotropic angular distributions of visible momenta, 
but differ in their m uD spectrum at higher y/s: the one of the WW fusion is always broad, while 
that of the Higgsstrahlung peaks strongly around m uD = mz as soon as y/s > m H + m z . In 
fact the two amplitudes interfere semiconstructively for y/s below and around the HZ threshold, 
but hardly so well above it. About 50% of the signal cross section around that threshold is 
associated with the WW fusion, this becoming about 10% at higher energies. 

The background consists of 4 classes of amplitudes, describing ZZ/7 cascade, WW fusion, 
W exchange and Z radiation processes. The ZZ and Z7 cascades dominate, the first leading 
beyond he ZZ threshold. Between the WW fusion and W exchange components there is a strong 
destructive interference and their coherent contribution is more than an order of magnitude 
smaller than the ZZ/7 one > while the Z radiation contribution is practically negligible. The 
overall behaviour over phase space is characterised by a peak around m vV = mz, two peaks in 
m b i due to the 7 — * bb and Z —» bb cascades, and a clear angular peak along the beam axis. 
The total background cross section (before cuts) lies between the mu = 90 GeV and m H = 100 
GeV signal cross sections throughout LEPII energies. The signal-background interference is 
negligible, even around mu = mz- 

Applying a set of minimal cuts prescribed by the LEPII groups, m b i > 50 GeV and m z — 25 
GeV < m vV < m z + 25 GeV, has a very different effect for y/s below and above the HZ 
threshold. At lower energies both the signal and the background cross sections are greatly 
reduced, although the first mostly due to the m uP cut while the second to the m bb cut. At 
higher energies both the signal and the background are reduced weakly (~ 10 %), the signal- 
background ratio improving only slightly. The effect of an angular acceptance of a b-tagging 
micro-vertex on the signal and the background were found roughly equal, rejecting about 10- 
15% of events for present type extended microvertices. The code exists also in the form of a 
MC generator. 
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Figure 1: Signal Feynman diagrams 
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Figure 2: Background Feynman diagrams 



Figure 3: Number of expected vvbb events from the background, the Higgs signal (m#= 90 GeV) 
and their sum, as a function of E C m for 1000 pb _1 . The signal+interference is also plotted, but 
is indistinguishable from the signal-only curve. The ratios of the signal-to-background and the 
total signal-to-Higgstrahlung contribution are also shown. QED ISR corrections are included 



Figure 4: Higgs signal vvbb cross section as a function of E CM for different HB masses and the 
background, without and with ISR corrections 



Figure 5: Total background vvbb cross section and subprocess contributions as functions of 
Ecm without and with ISR corrections 



Figure 6: Differential cross sections for the Higgs signal with mn 90 GeV (solid line) and for 
the background (dashed line) at Ecm=192 GeV. (ISR corrections included) 



Figure 7: Differential cross sections for the background at 192 GeV (solid line) and 175 GeV 
(dashed line) (ISR corrections included) 



Figure 8: Differential cross sections, at E C m=^^ GeV for the total background (solid line) and 
the subprocesses a) ZZ/7 production (dashed line) b) W fusion and exchange (dotted dashed 
line) and c) Z radiation (dotted line) (ISR corrections included) 



18 



Figure 9: Total background vvbb cross section and subprocess contributions as functions of 
Ecm, without and with ISR corrections, after cuts 



Figure 10: Higgs signal vvbb cross section as a function of Ecm for different Higgs masses and 
the background, without and with ISR corrections, after cuts 



Figure 11: The number of expected Higgs vvbb events for 1000 pb _1 as a function of E C m for 
m H =100 GeV, before and after cuts. The ratios of the signal-to-background and the total 
signal-to-Higgstrahlung contribution are also shown (ISR corrections are included) 



Figure 12: Percentage of events for Higgs (dashed and dotted dashed lines) and background 
(solid and dotted) outside the fiducial acceptance of a microvertex extending to cos# when 
demanding a) at least one taged jet or b) both jets tagged 
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